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ABSTRACT 

This theoretical study of the trajectories of charged particles 
supplements an experimental project of molecular-ionic rearmangements 
using the magnetic focusing properties of a beta-ray spectrometer, The 
experimentally measured magnetic field was analytically represented 
by a twelfth-order polynomial. This field is axially symmetric, but 
non-homogeneous otherwise. The particular particle of concern was 
the hydrogen ion, He. The trajectory of this particle was computed 
from a second-order differential equation, assuming values for the 
kinetic energy and initial angle of the particle, and the magnetizing 
current. The solution was obtained by numerical integration using a 
CDC 1604 digital computer. The distinctive feature of these calcu- 
lations, in contrast to those used normally for a beta-ray spectrometer, 


is the large scattering (initial) angle, about 45°. 
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1. Introduction 

This presentation is part of an experimental project presently being 
conducted at the U. S. Naval Postgraduate School to verify a classical 
theory developed by Bates, Cook, and Smith [1 ] of the capture of a 
light ion or atom from a target system by a fast projectile. 


Specifically, the theory involves a reaction of the type: 
H’'+CH, - H H+CH, 


and is essentially a sequence of binary collisions. The process 
ensues with the bombardment of methane gas by low energy protons, 
whereby the resulting rebound proton and escaping hydrogen ion join to 
form the product H, . The angle of departure (or scatter) is shown to 
be in the vicinity of 45 degrees, and is labeled @. 

Since the expected cross section will be quite low (<107° cm”), 
a detector could not be utilized, as is standard procedure, by placement 
at the scattering angle. [5 ] 

Therefore, the beta-ray spectrometer was chosen due to its prop- 
erty which will cause all ia ions of the same departure angle and 
energy to converge. This study determines that point of convergence, 
or focus. 

2. Discussion 

It is known from elementary electromegnetic theory that a magnetic 
force, a{v x B}, is produced on a charged particle moving at an angle 
to the direction of the field. This force causes the particle to sustain 
a change in velocity (direction only). The supposition is, for this 
field, that a particle which begins its motion on the axis of the field 
and at an angle with that field axis, will travel a curvilinear path and 
ultimately intersect the axis of the field for appropriate values of the 
field strength, energy and initial angle of the particle. This can be 
readily shown for a homogeneous field. 

The point at which the particle returns to the axis is called the 
intercept point, and the distance from the origin of the trajectory to 


the intercept point is the intercept distance (L). This is where the 


detector can be placed to analyze the particulate beam. The experimental 
usefulness is that if the intercept distance is known along with the 

field strength, and the energy of beam E,, then the initial angle may be 
evaluated, 

The theoretical problem, which is herein shown, is the evaluation 
of the intercept distance assuming values for the scattering angle. 

Thus, the following intermediate problems arose: 

A. A means to express the field, or more precisely, a continuous 
function which would represent the field for all points in the area of 
interest. 

B. A set of equations to specify the motion of the particle to 
incorporate the energy and charge of the particle, as well as the inter- 
action of the magnetic field. 

C. A numerical method to solve a differential equation which 
would result from the equations in (B) above. 

D. Interpretation and extrapolation from sets of trajectories to 


obtain these derivatives: 


(2 -- The change in intercept distance with 
ob/6, Ip respect to changes in energy 


25 
(2, Iz, Eg -- The change in intercept distance with 
respect to changes in angle 


ie 
(22) Eo, 8 -- The change in intercept distance with 
respect to changes in magnetization current 

3. The Magnet and its Field 

The magnet which will be used for this project was previously in 
the USNPGS beta-ray spectrometer. It is essentially composed of a 
normal section of a ‘right circular cylinder, and wound with cylindrically 
concentric layers of copper wire. The outer disks and cylindrical 
center section which form the container for the layers of wire are 
copper as are the cooling water tubes. The figure (Fig. 1) on the 


next page shows the magnet with its dimensions. 
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The experimental field strength of this magnet on-axis, H), appears 
as a Gaussian curve as shown in Fig. 2. The quantities in Appendix 1 
for the values of the field are shown in the positive z direction. The 
maximum intensity of the field on-axis is at the geometrical center of 
the magnet, at the point formed by the intersection of lines ii and jj. 

B has cylindrical symmetry about the line jj, and is symmetrical about 
the center plane. 

In order to more easily discuss the field, the system of cylindrical 
coordinates should be established at this point. The origin is defined 
to be at the intersection of the lines ii and jj. The line jj will be 
the z axis, and the radial distance, r, will be measured perpendicularly 
from it. The angular component, ¢, is as in (r, ¢, 2), a right hand 
coordinate system. 

Measurements were made on this field for garious values of field 
current by P. J. Kelly in May 1965 with an axial Hall Probe, and this 
data has been compiled in Appendix 1. The axial component of the 
field intensity was measured initially on the primary axis of the magnet 
(2 axfs) for a distance of 66 cm. on either side of the origin at 4 cm. 
intervals, This was done for various current settings between 2 and 
14 amperes, 

Also, the axial component of the field was measured off axis at 
four values of z, and for 8 values of ¢, for each value of z. The 
radius r was varied from zero to seven cm, 

The conclusions concerning the field were that: 

(1) The field strength is linear with current over the entire 
range of current settings within an accuracy of +0.1%. 

(2) The field is symmetrical on either side of the z = 0 
plane to +0.01%. 

(3) The field is azimuthally symmetric to +0.01%. 

To find the analytic expression for the field on axis, H,(z), where 
Zz assumes positive values only, a polynomial was fit to the data by 


the method of least-squares. This data is contained in Appendix 1. A 
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twelfth order polynomial was chosen, and provided sufficient precision 


in the form; 
H, (z) = B, + Bgz + Bgz* + Baz? + Boz* +... + Byyz? . (3-1) 


The thirteen coefficients were determined by the computer and a USNPGS 
subroutine LSQPOL, which is in PROGRAM CURVE included in Appendix 
6. 

At this point, a comparison was made of the experimental field 
values and those values which would result from the polynomial field 
expression. 

To make this eheck of the values of the field off-axis, the following 
steps were taken: 

(1) The vector potential was computed for all values of r and 2z 

using the derived polynomial. 

(2) The axial component of the field was derived, using the 

vector potential, and evaluated for various values of r and Z. 
This computation is shown in detail in Appendix 1. The agreement be- 
tween experimental and theoretical values was again quite good. 
4. Theory of the Trajectory 

The H," ion has a positive charge eqyal in magnitude to that of the 
electron, and mass essentially equal to twice the mass of a single 
proton. The mass of the single electron is ignored. 


The path of this particle can be represented by this equation: 


i 2_ ye ae aeal CA oA = - 
F Oa (os Sst a] +a [2 : aa 


where A is the vector potential, and k is defined as (mv/q). The 
prime above indicate differentiation with respect to z. A discussion 
of A, its derivatives, and k will follow. 

This equation for any charged particle can be derived from the 
equations of motion for that particle which begins its path on the 
primary axis of the magnetic field. This equation is derived in Appendix 


4. 
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The energy of the H, particle after scattering is E, and can be 
shown to be: 

E=2E, cos*® @ (4-2) 
where E, is the kinetic energy of the particle which is incident on the 
target methane. The angle @ is again the angle of scatter. 

A description of the collision process is as follows: 
The incoming proton strikes one of the hydrogens in the methane molecule. 
This is the first of the two binary collisions, and the only one which is 
assumed to be totally elastic. The initial particle moves off at an 
angle, called the departure angle; and the second particle (essentially 
a proton and an electron) moves towards the carbon atom at an angle 
equal to 90° - @. Since the carbon atom is assumed to be infinitely 
massive, compared to the hydrogen atom approaching it, this second 
particle is redirected in the second binary collision in the direction of 
the recoiling proton. 


A schematic of the collision process is as follows: 


He 
cS 
xy" 
2» 7 ‘) 
V eae 6g 
@)—»> — ee ee aes 
yy 
Pe 








carbon 
In the laboratory frame, the initial energy of the system is: 
Eo=$m, (¥,)° (4-3) 


and the final energy is: 


E = + (2m) (v")* (4-4) 
Conservation of momentum yields, for the first collision: 
P, =mv = mv' cos @ + mv" sin 6 (4-S) 
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and, 
Py =mv' sin @- mv" cos 6=0 (4-6) 


The last equation results in 

Vv" =V' tan @ (4-7) 
and substituting into the equation (455) 

> 


vi’ =vcos 6. (4-8) 


With this relationship between the velocities, the final energy after 


the two collisions becomes 


E = $ (2m) (v')* (4-4) 
= ¢ (2m) v® cos* @ 
= 2E, cos* @ 

b= 2E- eos 6 (4-2) 


If both collisions are elastic in nature, then it is quite easy to show 
that @ must be 45° . The y-component of momentum produces the 


equation 
es mv' sin @- mv' cos @= 0 (4-9) 
at the first collision. And it can only be satisfied when 


sin 6=cos @ 
Or, when me = 45": 
59. Expression for the Vector Potential 

The fact that the field is cylindrically symmetric carries with it 
the implication that there is no aZimuthal component, since the curl 
B=zero. This field can be described by a vector potential which is 
always in the azimuthal direction. 

Using the calculated expression for magnetic field on-axis, H,(z), 
the vector potential expansion 


is Pea peel AR eee 
A(r,2) ; H(z) 7 ey hell Uae ee Ha \2) 
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(S-1) 


is shown to be for any point (r,2): 


As = [G,]--(G,]+1G.]-— [G,.1+ ies 


384 18432 1474560 
ta [e414 | 
ee a eer ean ee 
176947200 ; P0727 29600 © 


n 
where (24 H,(z) is the nth partial derivative of the function 
Zz 


2 
H,(z) with respect to z, and where G, = H(z), and G, -(2} Halz); etc. 
Zz 


6. The Computer Solution 

It remains, at this point, to find the solution of the differential 
equation. The method chosen for this task was a process of numerical 
integration employing a digital computer. The subroutine INTEG-1 of the 
USNPGS Computer Library was selected because of its adaptability to 
this particular problem. The program in which INTEG-1 is contained 
is called PROGRAM ORBITS, and is listed as Program C in Appendix 6. 

The first four statements of this program (ORBITS) are standard 
according to the instructions for use of INTEG-1. The three comment 
cards that follow describe ENER, OKAY, and P, which are names given 
various quantities used in the solution. 

B(1) through B(13) are the coefficients used in the polynomial 
as obtained by PROGRAM CURVE. The functions F(1) through F(6), 
and G(l) through G(7) are explained in Appendix 5. 

VEC DZ is the symbol for the partial derivative of the vector 
potential with respect to z. VECDR represents the partial derivative 
of the vector potential with respect to r. VECPT is the vector 
potential. 

A second order differential equation can be reduced to a system 
of two first differential equations, Applying this technique to 


equation (4-1) yields 


oA oA oA oA 
Se } =_' [A oe, a. a [A pore (r*)* [A a a [A ae ] 
[ke = A] [ke - RP] [k? - At] Lie — Re 
(6-1) 
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eae eer (6-2) 
Letting 
x(1) =r and x(2) =r’ , 
the equations (6-1) and (6-2) become 
[x(2)]' = [x(2)] FUNC4 + [x(2)]® FUNC4 - [x(2)]* FUNC5 - FUNC5 
[x(1)]' = x(2) 
where 


of Sipe) 
FUNC4 — 


[k* - A*] 


FUNCS Sr 


ieee | 





Further, FUNC4 and FUNCS5 can be written in terms of FUNC1, 
FUNC2 and FUNC3 as follows: 


EUNCEL = ko =A" 
FUNC2 = A SAL 
Oz 
FUNC3 =A 9A 
or 
and 
FUNCI1 
FUNCS =_FUNC3 
FUNC1 


OKAY is obtained in the following way. This is the machine 
language expression for k* where the energy can be entered in eV. In 
the MKS system, k must have the unit of webers/meter. This is the 


a = ss 
same as A, since B has the units of webers/meter*, and B=7xA; 


k= Mv_ inunits of (kilogram)(meter/sec) or webers/meter. 


g (coulombs) 
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It would be desirable to have k expressed in terms of the energy 
of the particle, E, and in the units of eV, as a great savings in con- 
version would be realized. 

Using the non-relativistic expression for the energy E, the 


following is obtained: 


ME 
E =e “. p= (2mE)® 
m 


But, since k* is needed 


2 \ 
k? =P = oe = 2mb' if E‘ is expressed in eV. 
2 q 2 


And, k* becomes 
2(mass of H,") ” 
charge of electron 
Inserting the values of the constants, k* becomes the quantity OKAY: 


OKAY cae a ead } BE 
1.6021 x 10°” 


= (4.17576 x 10° JE" 


or 
OKAY = (4.17576E - 08)*(ENER), as expressed in the program. 

The method of solution employed by INTEG-1 is a fourth order 
Runge-Kutta numerical integration approximation of the ordinary differ- 
ential equations. 


The computation is performed in the MKS system of units. 
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TO USE PROGRAM ORBITS (the main trajectory program) 


This program is written in the language of FORTRAN-60. Its 
successful employment depends on the correct use of the data cards, 
which are the last eight cards of the program itself. These eight cards 
contain the input information, and the quantities desired as output. 

An example of the data cards may be seen on the last page of Program 
C in Appendix 6, 

Data card #1: This card is used for identification, and columns 2-32 
are available for this use. This title appears on all data output, as 
well as graphs. 

Data card #2: This card contains the number of "runs", or trajectories, 
to be processed with a maximum of 9. This number is entered in 
column 1 of the card. Note, however, that only one graph is possible 
for the entire program. 

Data card #3: The number of coupled first-order differential equations 
is specified by this card (maximum of 30). In this case, there are two 
and this number is posted, right justified, in columns 1 and 2, as 02. 
Data card #4: The initial and final values of the independent variable, 
z, over which the integration process is to be carried out. The inte- 
gration step size is also stipulated. The integration can be processed 


in up to 3 segments, each with a different step size, thus: 
zi—Sl, 7)-82,72 83.77, 


or 7] wg} 2 g7F ? 


or ZI—S ZF, 


where ZI refers to the initial value of z (in this case equal to zero), 
ZF refers to the final value of z, Zi refers to intermediate values of 
z, and Si are the step sizes. All of these numbers are entered in 
meters. The corresponding values, with decimal points, are placed 
in the above order in columns 1-10, 11-20, 21-30, etc. 

Data_card #5: Columns 1-10 contain the value of the energy E in 
electron-volts, again, with the decimal point. Columns 11-20 


contain the quantity P which adjusts the magnetic field according to 
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the base of 10 amperes magnetization current. For example, for a field 
of 7 amps, the quantity 0.70 would be entered anywhere in columns 
11-20. 

Data card #6: Here are the values of the initial conditions. The format 
is the same as for cards 4 and 5. The first value expressed is the 
initial value of r (equal to zero), and the next value is the slope of 


the function at that point, (os . This derivative may be considered 
dz 


to be also the tangent of the scattering angle, 6. 
Data card #7: This card controls data print-out. A blank card will 
suppress this output. Again, using the ten column sequence, this 
format is as follows. The first 8 columns contain the title for each 
variable output desired, and the last two columns contain the variable 
identifier. The variable identifiers are 00 for z, 01 forr, and 02 for 
r'. Anexample: Z DIST OOR DIST OI1SLOPE 02 
Data card #8: This last data card controls graph output. A blank 
card will suppress this output. Remember: that only ONE graph for 
each submitted program can be drawn and it corresponds to the first 
set of data. The information is contained in columns 1-20. Columns 
1-16 contain the graph title, and columns 17-20 indicate the variable 
identifiers. Columns 17-20 must always be 0100. The origin, 
direction of plot, and scale distances have already been set within 
the main program. 
Running Time: The time required for the computer solution is dependent 
on the step size and the total range of the independent variable. This 
running time has averaged about three minutes for a single trajectory. 
7. Results 

The outcome of this project was the creation of a single-package 
system which generates hydrogen ion trajectories in an axially 
symmetric field. 

Approximately 100 production trajectories were computed for 
values of energy E ranging from 97 eV to 450 eV; magnetizing current 


from 7 to 15 amperes; and scattering angle of 35-55 degrees. The 
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numerical output data and the graph plots showed these orbit con- 
figurations to be quite smoothly continuous and consistent. 

The error in evaluating the intercept distance L was due mainly 
in the uncertainty of reading the graph or interpolating the data print- 
out. These values of L ranged from 20 centimeters to about 80 centi- 
meters, with an estimated error of +1.0%. The maximum radial distance 
of the orbits was around 8 cm. 

Comparison of these computer trajectories with those obtained by 
the fine wire technique used by Kelly showed only fair agreement to 
within +10.0%. The disparities in the different sets of trajectories 
have not as yet been determined. 

This computer system was checked with a constant uniform 
(homogeneous) magnetic field against an analytic solution. The com- 
parative results were excellent (to within +0.01%). 

Based on a Standard Trajectory of E, = 300 eV, 

I, = Magnet current = 12 amperes, 


@ = 45 degrees: 


oL 

Oi =e 5 One d 

a cone cm/deg 
16 

a = -9,7 0.8 cm/amp 
B ie Eo 


3) 1, 8 0.20 +0.01 cm/eV 

The above derivatives show that as @and Ip are increased, Lis 
decreased. However, an increase in E, produces an increase in L. 
Each of the partial derivatives were obtained with respect to a change 
in a single parameter, with the other two parameters held constant. 
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APPENDIX 1] 
DATA ON THE MAGNETIC FIELD 

The magnet which provides this inhomogeneous field was formerly 
employed in the beta-ray spectrometer which was reconstructed here 
at the USNPGS in the period 1955-56. Although much was known about 
the magnet itself, there existed no data on the magnetic field until 
measurements were made by Kelly [5 Jin 1965. 

This data has been compiled and is herein available for the field 
strength both on- and off-axis; and also, for the theoretically computed 
as well as experimentally measured values in Gauss. 

In the diagram of the Magnetic Field Strength in this appendix, 
there are two values listed for each (r, z). The experimental values 
in the diagram are the upper ones listed in each square. 

The theoretically computed values are listed in squares and are 
the lower values. The method by which these values were computed 
is as follows: 

The field on-axis is given by: 


Ho(z) = B, + Boz + Bgz* + B,2° + Bez* + ... + Bygz” (3-1) 


(z expressed in meters) 
where: 


1) = 8.0498907E-02 
2)= 2.3351831E-02 
3) = -4,8664438E+00 
4) = 4.2887590E+01 
5) = -4.3626183E+02 
6) = 3.9271101E+03 
B( 7) = -2.2648265E+04 
B( 8) = 8.2393030E+04 
B( 9) = -1.9359723E+05 
B(10) = 2.9437706E+05 
B(11) = -2.8025472E+05 
B(12) = 1.5200646E+05 
B(13) = -3.5865075E+04 


B( 
B( 
B( 

B( 
B( 

B( 


Since A g is the only non-zero component of the vector potential, the 
axial component of the field strength, Hz , can be computed and compared 


to values obtained by the Hall Probe. 


Vas 


H , becomes ate (Cr Ag i} 


where Ag is equation (5-1). 


And, after performing the indicated operations, 





© ie 
Hy, a Malar Gz - Ge eres (11) 


fc. 
64 2304 


This is demonstrated in Appendix 6. 
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APPENDIX 2 
USEFUL FORMULAE IN CYLINDRICAL COORDINATES 


Conversion from X=r COS ¢ 
Rectangular y=r sin ¢ 
Coordinates Z=Z 


The Gradient 
7¢ - (28)? + 1 (ae 
or r\o¢ 


The Divergence 


alas ic ae ih ays a ere 


o> 
— 
Fi 
\e 
~ 
Tee 








The Curl 
TxA =l|0Az - rAd) | ¢ 
r] o¢ oz 
_| d0AZ _ B8Ar A 
or az 
, 1 [a@Ag) _ adr] 7 
r or o¢ 


The Laplacian 


7#@ = ifa_(; ae} 1{at.e). (326 
ee Urelg or lag? ora 
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APPENDIX 3 
POWER SERIES EXPANSION TO OBTAIN THE VECTOR POTENTIAL 
The well known expansion, equation (5-1) will now be derived. 
This expression describes the vector potential for an inhomogeneous 
magnetic field, which acts only in the azimuthal direction. The field 
on-axis, Hy (z), and its even derivatives can be computed. 
Basically, the derivation of this expansion is as follows: 


Since the field is constant and involves no currents in the region of 


the iield, 


And, 
B=7x &, 
so that the vector equation for which a solution is desired is: 
7x (Fx A)=0. 
Conditions for this field are: 
Bae = 03 
and A, = A, = 0, 


In cylindrical coordinates, 7 x B becomes 


a ee Be. 26a Th [Pe 28s nt [20Bg) 2B, 
i log az or az r or a¢ 


and each component must be equal to zero, 

















oBZ = dB 
or 3z (1) 
Br = i OAz] _ 8 (rA o ) rae _ o(rA A ) (2) 
r o¢ ee4 r OZ 
By = dAz dAr = 0 
Or 02 
Bz = I 0 (rA Pia ene 2 1 6 (rA , ) (3) 
8 or o¢ 1g or 
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From equations 1, 2 and 3, 


oO | 0A A SOP 20h} ek peed) ye 
or EE i ‘| a r a a| 0, 








and the partial derivative cae is zero as the coordinates are indepen- 
dent. 

Thus, 

oA bgone 2 a eee ey 

or? r or e a2 Bes 
Here, a solution is assumed in the form nz oe 4 a where 


the coefficients a, (z) are differentiable. To determine the coefficients, 








a = Dnap(z) ra-l ; 
and 

oe = Z n(n-1)ay (2) pn-2 
and 

34 = Za,‘ (z)rn 
and 

2-2 = La," (2) m 


Substituting in, and multiplying through by r* 


E(n)(n-1)ap(z)r" + Ona,(z)r® - Cap(z)r™ + La_"(z)rMt2 =0 


> {nfm-1) +n-1}ap(z)r® = - Da,"(2)mté 
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For n=0, the left hand side of the last equation reduces to -a,(z), and 


there exists no n=] term. Therefore, 


-ao (4) +2, {n(n-1) + n-l}ap(z)r™ + n= an" (z)rnt2 = 0 


-a,(2) + 3 {(n+2)(n+1)+n+1 } apnag (z)r +2 zy > a," (z)rn*e =, 
n=0 n=0 


-a,(z) + > ((n4t3) (nl) Pana (z)rn*2 4 3 a,"(z)jrnté - 0 
n=0 n=0 


-a(z) + 2B { ( (n+3)(nt+1) Jaj4o (2) +a,"(z)}r Me = 0 
n=0 


The following must then be true, as for all power series expansions: 


ao(z) = 0, hence 
(n+3)(n+l)agig (Z) + ay"(z) = 0 
Siena) ee ee 

nee” ar) 


Whence, a recursion relationship yields the important needed para- 


meters. Obviously, there can be no terms for n even. 


For n= 1 ag(z) = 2a." (2) 
n= "3 a,(z) = —zaa (2 , ete. 
24 
As Alr-z) = > =e 4 
n=0 


=a,(z)r + ag(z)r? + ag(z)r° + aj(z)r? +... 
Making the substitutions, 
tt i i 1 ore 
A =ayle)r-ay"(2)= + (4, J) ay" "ah =... 
And, to evaluate a,(z), solve for Bz using 7 x A, 


6 (rA 
or 


2 





1 4 ia) 6 tboue 
Bz =7° = 2a; (2) - 5 ay (z)r* + 54g A (z)r* - 2... 
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Knowing that B-_ is equal to H,(z) for r = 0, hence, 
H(z) = Bz (r=0) = 2a, (2), az) SrtoH.(z) 


and the expression for the Vector Potential is obtained as 


r c oa 
A(r,z) = 7a Hp(2) ai fle 2) + 364 Hee ey = ete 
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APPENDIX 4 
THE DIFFERENTIAL EQUATION DEVELOPMENT 

The essentials of this operation involve simply applying Newton's 
second law in cylindrical coordinates to the magnetic force acting on 
the particle. Finally, time and the azimuthal variables can be elim- 
inated by means of the two motion constants, energy and the canonically 
conjugate angular momentum. Since the particle will begin its tra- 
jectory on the z-axis, the canonically conjugate angular momentum is, 
in fact, zero. 


Applying Newton's 2nd law, 
Pag = q(v xB) = m(a) 
where 


In cylindrical coordinates, the velocity vector is 
= ea 4 A 7 A 
Vv=errt+rg@qgt+ ZZ, 
and, the acceleration vector is 
a= (§-1 o?)T + (2r g + ré)¢ ee Ale? 
Knowing the vector potential, the magnetic field B and the 


quantity Vv Xx B can be found. 


ee Se ee eee nee * ae ee ee 5 
Beyx R= - Ena Zaqi % 2-(a2| sz (At 3 


= A 
But; Kh = Ag @; 


a 


Hence: 


S INe 


<} 
x 
wf 
ll 
| e 
ole 
Le J 
> 
oa 
| eee | 
n> 
~ 
! 
oa 1 (hae 


Hd) 
+> Bade. 


Equating each orthogonal cylindrical component, the following is 


3 A 
24 rA | g 


obtained, 
m@f-1$?) = a(d S—rA6), (4) 
m(2rg +r) =4 |-é Z(ra d) - s sean ") : (S) 
" . od 
m(z) =q(¢ re o). (6) 


These three equations can be reduced further to: 


* Po) 
+ Q@ ee o), 


and 


eo ath ior! 
“ar m2) = ar @ a AO). 


From the middle equation above, it can be shown that the quantity 
conserved is the canonically conjugate angular momentum about the 


symmetry axis. Thus, 





Pm = mre gt arg, 
4 = Po - QrAg@ 
2 
mr 
and, 
pS Ga eee: Pe Pc oA g 
dt om) m ar ane gr* or 


Replacing the independent variable of time (t) with the distance on 


the z-axis (z) yields 


a [> F *\2 2 Po -qrAg |* 
v* = [ =| + 7 (2)- + fr ee oe (7) 
and 2 : 
2 g ie 2 
(z)? = \ — me [ 2: 2 | 





ds . as 
arm) = ile a) i 
a eo A ° A o A 
Vv=eHarrtrg@g@tzz , 
v= yey = (1) sea) 2 2) 
and 
pe @ ; 
dz ; 
$= Po - aA 
mr* 
Hence: 
ar \ = ‘ poe orn. <= 
2. |({at 2 Cc 
V es + 1 (z)? + r® eee | 
_ ar \° *\5 q* Pc 2 
= } ts) (z)* + =e aan 
Further: 


ce . 2 = {2} [ee -a}? 
(z)* = m qr (3) 


35 


anda: 


denise dr i cies me ae Po _ Pc oA 
mdz? (z)" + dz E dz a| Sean {2 : q + ae (9) 


o@ [| = = {2a} e 2 sf (10) 


Combining the remaining equations 8, 9, and 10, and then allowing 


that 








Po = 0 
and 
ae: pe ; 
gq 


the following equation results: 


7! F . aA ¢ aAg |_ 
(4-1) 


where the dots throughout this discussion indicate derivatives with 
respect to thme, and the primes indicate derivatives with respect to 


the variable z. [4 ] 
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APPENDIX 5 
DEVELOPMENT OF THE VECTOR POTENTIAL 
To adequately describe the field at every point, and for use in 
the general trajectory previously derived, equation (4-1), the following 


quantities are needed: 


B = The Vector Potential, 
3A | | wa | 
a = The first partial derivative of A with respect to Zz, 
and 
3A | | a | 
= = The first partial derivative of A with respect to r. 


The vector potential, found by the power series expansion in Appendix 


a me RS: 
i i Ee reg 
A a 2 [G,]- 16 [Ge] = 384 [G3] = 18432 [G,] 5 ee 
= ‘ais 
29727129600 [Gr] , 


where, for the purposes of programming, the G functions are defined 
as follows: 
G, = He(z) 
_ 9° H(z) 
G2 as 3 Zz? 
4 


a 
3 3 z* 


aoe o° H, 
5 aoe 


= 9 Hy 
G, = a zie 


o7 


The first partial of A with respect to r becomes, 


= 





- a 9r® 
= #[G, 1-4 —[G,]+ 2 a 7 (Ga1- er [Gal+ Ta74560 bGs} 
Li 3c. 
- 176047200 ts] + 39727129600 LGr! 


The first partial of A with respect to z is a little more tedious, but 


can be shown to be, 


3A 
OZ 








= +Ih1]- ie + san a) ee ae Gg 


where the F functions are defined as follows (again for programming 


convenience), 


ee OG, 
pla gear dz 
aG 
c= 
: az 
_ Ge, 
Peja eee 
3 oe 
8G 
a 
; OZ 
0G. 
|e r= Eee 
5 52 
Fe = 0G. 
OZ 
F, = 0 
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APPENDIX 6 
PERTINENT MACHINE PROGRAMS 
The pertinent computer programs used in the solution of this 


problem are listed in this appendix in the following manner: 


Program A - Least-squares polynomial generation by the PROGRAM 
CURVE, and contains 7 pages. 

Program B - Theoretical off-axis field check labeled PROGRAM 
WATCHOUT, and contains 2 pages. 

Program C - The main trajectory program labeled PROGRAM ORBITS, 


and contains 15 pages. 
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APPENDIX 7 
GRAPHS OF VARIOUS TRAJECTORIES 


This appendix contains these four figures which show trajectories 


as produced by PROGRAM ORBITS. 


Figure 3. 
Figure 4, 
Figure 5. 
Figure 6, 


The Standard Trajectory 
Variation in Energy 
Variation in Magnetic Field 


Variation in Initial(scattering) Angle 
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